Next Generation Relativistic Models 



R. J. Furnstahl 

Physics Department, The Ohio State University, Columbus, OH 43210 
furnstahl . l@osu . edu 



1 Introduction 

The current generation of covariant mean-field models has had many successes 
in calculations of bulk observables for medium to heavy nuclei, but there 
remain many open questions ^E]- New challenges are confronted when trying 
to systematically extend these models to reliably address nuclear structure 
physics away from the line of stability. In this lecture, we discuss a framework 
for the next generation of relativistic models that can address these questions 
and challenges. We interpret nuclear mean- field approaches as approximate 
implementations of Kohn-Sham density functional theory (DFT), which is 
widely used in condensed matter and quantum chemistry applications 
We look to effective field theory (EFT) for a systematic approach to low- 
energy nuclear physics that can provide the framework for nuclear DFT [S]. 

We start with the key principle underlying any effective low-energy model 
and then describe how EFT's exploit it to systematically remove model depen- 
dence from calculations of low-energy observables. Chiral EFT for few nucleon 
systems is maturing rapidly and serves as a prototype for nuclear structure 
EFT. The immediate question is: why consider a relativistic EFT for nuclei? 
We show how the EFT interpretation of the many-body problem clarifies 
the role of the "Dirac sea" in relativistic mean-field calculations of ground 
states and collective excited states (relativistic RPA). Next, the strengths of 
EFT-motivated power counting for nuclear energy functionals is shown by 
the analysis of neutron skins in relativistic models, which also reveals weak- 
nesses in current functionals. Finally, we propose a formalism for constructing 
improved covariant functionals based on EFT/DFT. We illustrate the basic 
ideas using recent work on deriving systematic Kohn-Sham functionals for 
cold atomic gases, which exhibit the sort of power counting and order-by- 
order improvement in the calculation of observables that we seek for nuclear 
functionals. 
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2 Low-Energy Effective Theories of QCD 

If a system is probed with wavelengths small compared to the size of char- 
acteristic sub-structure, then details of that substructure are resolved and 
must be included explicitly (e.g., see Fig.^for a nucleus). On the other hand, 
a general principle of any effective low-energy theory is that if a system is 
probed or interacts at low energies, resolution is also low, and what happens 
at short distances or in high-energy intermediate states is not resolved [HI [7] . 
In this case, it is easier and more efficient to use low-energy degrees of free- 
dom for low-energy processes (see Fig. EJ. The short-distance structure can 
be replaced by something simpler (and wrong at short distances!) without 
distorting low-energy observables. This principle is implicit in conventional 
nonrelativistic nuclear phenomenology with cut-off nucleon-nucleon (NN) po- 
tentials (although rarely acknowledged!). There are many ways to replace the 
structure; an illuminating way is to lower a cutoff A on intermediate states. 
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Fig. 3. Tlie equation for the T-matrix witli cut-off potential Va and replacement of 
the effects of high q intermediate states with a contact interaction. 
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Fig. 4. The ^So Va &t A = 2.1 fm"^ (symbols) for many xVdof « 1 potentials 0. 
2.1 Renormalization Group and NN Potential 

We illustrate the general principle by considering nucleon-nucleon scattering 
in the center-of-mass frame (see Fig. The Lippmann-Schwinger equation 
iterates a potential, which we can take as one of the x^/dof « 1 potentials (as 
in Fig. 21). Intermediate states with relative momenta as high as q = 20 fm"'^ 
are needed for convergence of the sum in the second term. Yet the elastic 
scattering data and the reliable long-distance physics (pion exchange) only 
constrain the potential for q < 3fm~^ [5]. 

We can cut off the intermediate states at successively lower A; with each 
step we have to change the potential Va to maintain the same phase shifts. 
This determines the renormalization group (RG) equation for Va |H|. We see in 
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Fig.0]that at A — 2.1 fm^^, the potentials have aU collapsed to the same low- 
momentum potential ("Viowfc")- We emphasize that this potential still repro- 
duces all of the phase shifts described by the original potentials. The collapse 
actually occurs for A between about 500 and 600 MeV (or 2.5-3 fm~^), which 
implies model dependences at shorter distances. Note that the collapse does 
not mean that the short-distance physics is unimportant in the S channels; 
rather, only the coarse features are relevant and a hard core is not needed. 
(For much more discussion of Mowfc, including plots of the potentials and 
phase shifts in all relevant channels, see [5].) 

The shift of each bare potential to T^owfc is largely constant at low mo- 
menta, which means it is well represented just by contact terms and a deriva- 
tive expansion (i.e., a power series in momentum) This observation illus- 
trates explicitly that the short-distance physics can be absorbed into local 
terms; this is the essence of renormalization! Further, it motivates the use of 
a local Lagrangian approach. In an EFT, by varying the cutoff (or equiv- 
alent regularization parameter), we shift contributions between loops and 
low-energy constants (LEC's), just like the shift between the high-lying in- 
termediate state sum and the potential Va- The long-range physics is treated 
explicitly (e.g., pion exchange) and short-distance interactions are replaced 
by LEC's multiplying contact terms (including derivatives). 

2.2 Effective Field Theory Ingredients 

The low-energy data is insensitive to details of short-distance physics, so we 
can replace the latter with something simpler without distorting the low- 
energy physics. Effective field theory (EFT) is a local Lagrangian, model- 
independent approach to this program. Complete sets of operators at each 
order, determined by a well-defined power counting, lead to a systematic ex- 
pansion, which means we can make error estimates. Underlying symmetries 
are incorporated and we generate currents for external probes consistent with 
the interactions. The natural hierarchy of scales is a source of expansion pa- 
rameters. 

The EFT program is realized as described in [7|: 

1. Use the most general Lagrangian with low-energy degrees of freedom con- 
sistent with global and local symmetries of the underlying theory. For few- 
nucleon chiral EFT (our example here) , this is a sum of Lagrangians for 
the zero, one, and two-l- nucleon sectors of the strong interaction: 

>Coft = Cttt, + C-nN + CnN ■ (1) 

Chiral symmetry constrains the form of the long-distance pion physics, 
which allows a systematic organization of pion terms 0. 

2. Declare a regularization and renormalization scheme. For the few nucleon 
problem, the most successful approach has been to introduce a smooth 
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cutoff in momentum. The success of tlie renormalization procedure is re- 
flected in the sensitivity to the value of the cutoff. The change in ob- 
servables with a reasonable variation of the cutoff gives an estimate of 
truncation errors for the EFT |1(J| . 
3. Establish a well-defined power counting based on small expansion param- 
eters. The source of expansion parameters in an EFT is usually a ratio 
of scales in the problem. For the strong interaction case, we have the mo- 
menta of nucleons and pions or the pion mass compared to a characteristic 
chiral symmetry breaking scale , which is roughly 1 GeV (600 MeV is 
probably closer to the value in practice). The nuclear problem is com- 
plicated relative to tttt or ttN by an additional small scale, the deuteron 
binding energy, which precludes a perturbative expansion in diagrams. 
Weinberg proposed to power count in the potential and then solve the 
Schrodinger equation, which provides a nonperturbative summation that 
deals with the small scale. The NN potential is expanded as: 

oo 

with Q a generic momentum or the pion mass, and (see [T] for details) 

A + 2(L-C)+^y,(d, + /,/2-2) , (3) 

i 

where the topology of the corresponding Feynman diagram determines v. 

The organization of Vnn hi powers of Q is illustrated in this table (solid 
lines are nucleons, dotted lines are pions): 
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The table indicates what gets added at leading order (LO), next-to- leading 
order (NLO), and next-to-next-to-leading order (NNLO or N^LO). The low- 
energy constants at LO and NLO, which are the coefhcients of the AN contact 
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terms (nine total), are determined by matching to the phase-shift data at 
energies up to 100 MeV; phases at higher energies are predictions. (There is 
also input from vrA^ scattering through C^n at NNLO.) 
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Fig. 5. Representative phase shifts at LO (dots), NLO (dashes), and NNLO (solid) 
as a function of lab energy (in GeV) |1U|. 



In Fig. [SJ the systematic improvement from LO to NNLO is evident |10| . 
This type of systematic improvement is what we seek for the many-body prob- 
lem. Recent calculations at N'^LO by Entem and Machleidt pd, and by Epel- 
baum et al. |12| show continued improvement. The work by Epelbaum et al. is 
particularly noteworthy in providing error bands based on varying the cutoff. 
All of the phase shifts are consistently predicted within the estimated trun- 
cation error. Entem and Machleidt fine-tune their potential to achieve x'^ ~ 1 
for the phase-shift data up to 300 MeV lab energy. This puts the potential 
on the same footing as conventional potentials, but sacrifices the controlled 
systematics of the EFT. 

An important feature of the chiral EFT is that it exhibits naturalness. 
That means that when the relevant dimensional scales for any given term are 
identified and factored out, the remaining dimensionless parameters are of 
order unity. If this were not the case, then a systematic hierarchy would be 
in jeopardy. The appropriate scheme for low-energy QCD, which Georgi and 
Manohar ^] called "naive dimensional analysis" (or NDA), assigns powers 
of /tt = 93 MeV and to a generic term in the Lagrangian according to: 

C,..-c,m.[-j^) [j-J [-^) f.A^, (4) 

where m, and n are integers. 

The dimensionless constants resulting when the LO and NLO constants 
are scaled this way |1()| are given in the table below for cutoffs ranging from 
500 to 600 MeV. We see that 1/3 < Q„i„ ^ 3 in all cases, with one exception, 
which implies they are natural (and the expansion is under control). The one 
exception is f^Cr, which is unnaturally small. This is often a signal that there 
is a symmetry, and in this case a corresponding symmetry has indeed been 
identified: the Wigner SU{4) spin-isospin symmetry |1(J| . 



Next Generation Relativistic Models 7 



Cs 


-1.079... - 0.953 


Ct 


0.002 . . . 0.040 
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0.403... 0.281 
2.846 . . . 3.410 
-1.929... - 1.681 
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2. 029... 2. 251 
-0.364 ... - 0.428 
-0.728 ... - 0.668 



We can also explore the connection between the chiral EFT and con- 
ventional potentials by assuming that the boson-exchange potentials provide 
models of the short-distance physics that is unresolved in chiral EFT (except 
for the pion). Thus this physics should be encoded in coefficients of the contact 
terms. We can reveal the physics by expanding the boson propagators: 

^ V + Xv^ + 
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By comparing to Q with m ~ A^, we see that g ~ A^/f^^, which explains 
the large couplings for phenomenological one-boson exchange! 
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Fig. 6. Comparison of EFT and one-boson-exchange coefficients |10|. 



In Fig. IHl coefficients from the chiral EFT are compared to those obtained 
from such an expansion applied to several boson-exchange potentials |1U| . The 
semi-quantitative agreement is remarkable and shows that there are strong 
similarities between the chiral potential and standard phenomenology. It is 
not yet clear, however, whether one can conclude that the potentials contain 
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reasonable models of short-distance physics; the agreement may just reflect 
the fact that they all fit the same data. 

In summary, the chiral EFT for two-nucleon physics (and few-body nuclei) 
is maturing rapidly. The systematic improvement and model independent na- 
ture is compelling. We seek the same characteristics for our description of 
heavier nuclei using covariant energy functionals. 

3 Relativistic versus Nonrelativistic EFT for Nuclei 

Why should one consider a relativistic effective field theory for nuclei? Let us 
first play devil's advocate and argue the contrary case. The relevant degrees 
of freedom for low-energy QCD are pions and nucleons (at very low energy 
even pions are unresolved). Nuclei are clearly nonrelativistic, since the Fermi 
momentum k-p is small compared to the nucleon mass. The nonrelativistic NN 
EFT described above is close to the successful nonrelativistic potential and 
shows no (obvious) signs of problems. 

In the past, a common argument against covariant treatments was that 
they intrinsically relied on "Z graphs," which implied NN contributions that 
were far off shell. The claim was that those should really come with form 
factors that would strongly suppress their contribution I14| . This argument 
does not hold water in light of the effective theory principle we have high- 
lighted: the high-energy, off-shell intermediate states may be incorrect, but 
that is fine, since the physics can be corrected by local counterterms. (Note: 
this only works if we include the appropriate local operators!) The modern 
EFT argument against Z graphs is different: since they are short-distance de- 
grees of freedom, they should be integrated out. Indeed, it is often said that 
the heavy meson fields (tr, w, and so on) should be integrated out as well. 

This argument is sometimes advanced as a matter of principle, but the 
underlying reason is the need for a well-defined power counting. The early 
experience with chiral perturbation theory using relativistic nucleons was that 
chiral power counting is spoiled by unwanted factors of the nucleon mass that 
come essentially from Z graphs 1151 . Integrating out heavy degrees of freedom 
moves these large scales to the denominators and then all is well |16| . 

We'll resolve the issue of power counting in the next subsections. But first 
a brief (in the legal sense) on the side of the (covariant) angels. Arguments 
in favor of covariant approaches to nuclear structure are collected in |S] and 
|17|: this is merely a skeletal summary. For nuclear structure applications, 
the relevance of relativity is not the need for relativistic kinematics but that 
a covariant formulation maintains the distinction between scalars and vec- 
tors (more precisely, the zeroth component of Lorentz four- vectors). There is 
compelling evidence that representations with large scalar and vector fields 
in nuclei, of order a few hundred MeV, provide simpler and more efficient 
descriptions than nonrelativistic approaches that hide these scales. The domi- 
nant evidence is the spin-orbit splittings. Other evidence includes the density 
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dependence of optical potentials, the observation of approximate pseudospin 
symmetry, correlated two-pion exchange strength, QCD sum rules, and more. 

3.1 Historical Perspective: Relativistic Hartree Approximation 

From the first applications of relativistic field theory to nuclear structure 
("quantum hadrodynamics" or QHD), the Dirac sea had to be considered. 
Indeed, it was often hailed as being new physics missing from the nonrela- 
tivistic description. The simplest (and seemingly unavoidable) consequence 
of the Dirac sea was an energy shift from filled negative-energy states in the 
presence of a scalar field 4> that shifts the effective nucleon mass to ill* (i.e., 
a form of Casimir energy, see Fig. for a schematic of shifted nucleon poles): 

(5i/ = -^[(k^ + i\r2)i/2 _ _^m2)1/2] ; M*=M-g,(f, (5) 

kA 

This sum is divergent. Two paths were taken to deal with the divergence. 
The "no-sea approximation" simply discards the negative-energy contribu- 
tions (and therefore SH) with a casual argument about effective theories (we 
fix up this argument in EFT language below). 
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Fig. 7. Pole structure for the single-particle propagator. 



The other approach was to insist on renormalizability of the QHD theory. 
The usual prescription eliminated powers of (j> up to cj)'^ in SH, leaving a finite 
shift in the energy in the "relativistic Hartree approximation" (RHA) of ^Hl 

1 M* 
5!^5 ,^if , 6!k6 ,,if 

In accordance with Georgi-Manohar NDA, we have in the second line scaled 
(f) with a factor of /^r and introduced combinatoric factors and dimensionless 
couplings K5 and kq of order unity (which absorb factors oi gsfn I M k, 1). The 
overall factor of M"*, however, is much larger than the corresponding factor of 
f^A^ in NDA. This is a signature that the power counting will not be correct. 

In practice, the (jr' and higher terms are large and drive the self-consistent 
effective mass M* too close to M, which results in too small spin-orbit split- 
tings (the raison d'etre of the relativistic approach!). Furthermore, the loop 
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Fig. 8. Contributions to tiie energy per particle from terms in the RHA energy 
functional proportional to (j>" ^M- Two values of M* /M are considered (diamonds 
and squares). The X's and error bars reflect natural NDA estimates. 



expansion is a phenomenological disaster (two-loop corrections qualitatively 
change the physics and there is no sign of convergence [201 )• The factor in 
(O also badly violates Nc counting that we would expect in a low-energy the- 
ory of QCD. {Nc is the number of colors.) With associated with a meson 

mass, oc (iVj" while A^f^ oc {l)^{Nc^^)'^ cx {N^y [H]. 

In general, the RHA energy functional exhibits unnaturally large contri- 
butions to 0" terms, as seen in Fig. |H1 A natural effective theory would have 
decreasing contributions, so that 0^ and higher terms would have only mi- 
nor effects (easily absorbed elsewhere). With a value of M* /M « 0.6, which 
is required for realistic spin-orbit splittings, the energy/particle from the (f)^ 
term is two orders of magnitude too large. The system adjusts M* upward to 
reduce this contribution, but the consequence is small spin-orbit splittings. 

3.2 Pov^rer Counting Lost / Power Counting Regained 

Gasser, Sainio, and Svarc first adapted chiral perturbation theory (ChPT) 
to pion-nucleon physics using relativistic nucleons |15| . In contrast to the 
pion-only sector, however, the loop and momentum expansions did not agree 
and systematic power counting was lost. The heavy-baryon formulation was 
introduced to restore power counting through an expansion in l/M, which 
means a nonrelativistic formulation jl6l. 

In 1996, Hua-Bin Tang wrote a seminal paper |23 in which he said 
". . . eft's permit useful low-energy expansions only if we absorb all of the 
hard-momentum effects into the parameters of the Lagrangian." The key ob- 
servation he made was that: "When we include the nucleons relativistically. 
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the anti-nucleon contributions are also hard-momentum effects." Once the 
"hard" part of a diagram is absorbed into parameters, the remaining "soft" 
part satisfies chiral power counting. 

There are several prescriptions on the market for carrying out this pro- 
gram. Tang's original prescription (developed with Ellis) involved the expan- 
sion and resummation of propagators |22I ■ This basic idea was system- 
atized for ttN by Becher and Leutwyler [2'A\ under the name "infrared regu- 
larization." More recently, Fuchs et al. described an alternative prescription 
using additional finite subtractions (beyond minimal subtraction) in dimen- 
sional regularization j^^. Goity et al. [23 and Lehmann and Prezeau 
extended infrared regularization to multiple heavy particles. The result for 
particle-particle loops reduces to the usual nonrelativistic result for small mo- 
menta while particle-hole loops in free space vanish identically. 

3.3 Effective Action and the No-Sea Approximation 

At the one-loop level, we will be able to apply the necessary subtractions 
without specifying the regularization and renormalization prescription in de- 
tail. It is convenient to use an effective action formalism to carry out the 
EFT at finite density. The effective action is obtained by functional Legendre 
transformation of a path integral generating functional with respect to source 
terms. This is analogous to Legendre transformations in thermodynamics. We 
assume some familiarity with effective actions; see [23 EH EH EOl to learn 
more. 

QHD models with heavy-meson fields correspond to effective actions with 
auxiliary fields, introduced to reduce the fcrmion integration to just a gaussian 
integral. We denote the effective action r[(j>,V'^], with classical fields and 
y (suppressing all other fields) (20j. An alternative is to work with a point 
coupling model, in which case the effective action will be a functional of the 
scalar nucleon density and the vector current. 

There will always be a gaussian fermion integral, which yields a deter- 
minant that appears in the effective action as the trace over space-time and 
internal variables of the logarithm of a differential operator: 



In the point-coupling case, optical potentials that are functions of the density 
and current appear directly rather than through meson fields, but the form is 
the same as ((HJ. This Tr In depends on /i, the nucleon chemical potential. 
The derivative expansion of the Tr In at /i = takes the form : 



Tr ln{i^ + ^7° - M* ~ g^J) = Tr InG^^ . 



(8) 




+ -Z2,(0)(n0)2 + ...] , 



(9) 
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which shows that this is a purely local potential in the meson fields. This 
means it can be entirely absorbed into local terms in the Lagrangian. That is, 
by adjusting the constants appropriately, this contribution to the energy will 
not appear. We specify a consistent subtraction at a specific /i, which means 
removing a set of local terms (that are implicitly absorbed by parameter 
redefinition). The obvious choice is /i = 0. 

We can accomplish this in practice at /i 7^ simply by performing a 
subtraction of the Tr In evaluated at /i = 0. Note that this is not a vacuum 
subtraction, because it depends on the finite-density background fields 
and V^{x). It is simply a choice for shifting "hard" Dirac sea physics into 
the coefficients. We emphasize that the same coefficients in the derivative 
expansion are subtracted for any background fields. For the ground state, 
this will simply remove all explicit evidence of the negative-energy states. 
But there are fixed consequences for treating linear response (RPA) that still 
involve negative-energy states 051 IMI 1^ . 

For the ground state, field equations for the meson potentials are ob- 
tained by extremizing r[(l),V^^]. These equations determine static potentials 
00 (x) and Vb(x) for the ground state; the corresponding Green's function 
is the Hartree propagator, denoted Gh- Then /^[0o(x), Vb(x)] is propor- 
tional to the zero-temperature thermodynamic potential Q = E ~ iiN. Both 
+iTT \nGJl^{^l) and -FiTr InG^i(O) are diagonal in the same single-particle 
basis ilJa{^)e^'^^° , where the '!/'a(x)'s are solutions with eigenvalues to the 
Dirac equation in the static potentials. Thus the subtraction is simple (T is a 
constant time that cancels out): 

r/T = -n= [-iTr InG^i(^) -^iTr InG^i(O) + • • •] /T 

= ^(A*-e„)- ^(-6.) + ---= (a*-^-)' (10) 

where we have omitted the parts of F (and fl) depending only on meson 
potentials, and the vacuum subtraction for the baryon number |35l I35j . Thus 
the contribution to the energy is simply a sum over occupied, positive-energy, 
single-particle eigenvalues. This is the "no-sea" prescription! 

Similarly, we can calculate the scalar density (the expectation value of -0?/;) 
from the effective action, which yields 

p,(x) cx 6rl6{gs4>{^)) =^ Tr Gh(m) - Tr Gh(0) . (11) 

The subtraction removes the contribution from the sum over negative-energy 
states. This is indicated diagrammatically as : 
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where "/i" denotes holes and "— " denotes negative-energy states. Thus, we 
recover the "no-sea approximation" for the ground state. What is the corre- 
sponding resuh for RPA excitations? 

Consider /^[(/f, y^] evaluated with time- dependent fluctuations about the 
static ground-state potentials: 

= 00 (x) + ^{x) , V^{x) = yo(x) (5^0 + V^'ix) . (12) 

We can access the linear response directly from the effective action, because F 
is the generator of one-particle-irreducible Green's functions [22! ■ The linear 
response is dictated by the second-order terms in powers of (/>, 1^^. So we 
expand 

r = -iTT lnG"^(^) + iTr lnG"^(0) + ... (13) 
in (this generates the "polarization insertion"): 

HGh' + 9s4>) = HG-H% • • - ^Gff {gs4>) Gh (.9.0) + ...]• (14) 
The two Gij's mean that we generate rings with Hartree propagators: 

~iT^\uG-\^x)^ pl^h + pI^- -iTrlnG-i(O) =^ +1^- (15) 
but different spectral content. Combining these contributions as in (fU^ : 







We see that the consistent prescription for the RPA is to include particle-hole 
pairs {ph) and pairs formed from holes and negative-energy states (h—). 

The lesson from effective field theory is that we must absorb hard- 
momentum Dirac-sea physics into the parameters of the energy functional. 
This explains the successful relativistic mean-field and RPA "no-sea" phe- 
nomenology. The QCD vacuum effects are automatically encoded in the fit 
parameters; there is no need to refit at different densities. In principle, to 
absorb all of the short-distance physics we need all possible counterterms. 
In practice, we rely on naturalness to make truncations (which justifies the 
success of conventional models with limited terms). 

The effective action formalism is a natural framework for systematically 
extending the one- loop treatment. By embedding the "mean- field" models 
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within the density functional theory (DFT) framework, we reahze that correla- 
tions are approximately included with the present-day mean-field functionals, 
which can be generalized to systematically improve them according to EFT 
power counting. Note that the effective action formulation leads naturally to 
time-dependent DFT as well. To carry out this program, a generalized infrared 
regularization scheme (or equivalent, possibly with a cutoff) is needed. 



4 Toward More Systematic Energy Functionals 

The energy functionals of relativistic mean-field models need improvement in 
a number of areas. As stressed in our first discussion of EFT's, long-range 
physics must be included explicitly. However, chiral physics as manifested in 
long-range pion interactions is not included explicitly in mean-field models and 
the treatment of long-range correlations is unlikely to be adequate. Pairing is 
still more an art than a science. Isovector physics in general has problems, as 
asymmetric (iV ^ Z) nuclear matter is poorly constrained. Here we focus on 
the last of these. 

EFT- inspired power counting (in the form of NDA [37]) can be applied to 
mean-field energy functionals to quantify the isovector constraints. Estimates 
of individual contributions in a given functional to the energy per particle of 
a finite nucleus can be made using NDA rules to associate powers of f-^ and 

with scalar (ps), baryon (pu), isovector (pa), and tensor (s) densities and 
gradients: 

^ ^ (Pb) ^ (VPb) ^ (3) ^ Z-N (fa) 
Ps,P.^^ Vfa-^ P^^^JT^ (16) 

and then using local density approximations to estimate the densities. Typical 
results are shown in Fig. El The hierarchy of isoscalar contributions is clear 
and shows that three orders of parameters can be constrained (contributions 
below about 1 MeV are not constrained). In contrast, the isovector parameters 
are weakly constrained by the energy; essentially only one parameter is deter- 
mined. Naturalness implies that the (Vps)^ contribution to observables such 
as the skin thickness r„ — rp is small (see |38|'). so we can focus on terms in 
the functional that contribute for uniform systems. This leads us to examine 
the symmetry energy in uniform matter. 

4.1 Power Counting and the Symmetry Energy 

The relationship between the neutron skin in ^'^'^Pb and the symmetry energy 
in uniform matter has been studied in detail [351 HOI IS HHI BH- We can 
frame the discussion in terms of the energy per particle E{p, a) in asymmetric 
matter, which can be expanded as: 

N — Z 

E{p, a) = E{p, 0) + S2{p)a^ + S4{p)a^ + ■■■ a = (17) 
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Fig. 9. Estimates in finite nuclei using averages of densities and results from two 
models (at left and right of error bars) |17| . 



where p is the baryon density and a is an asymmetry parameter. The a = 
term has the usual expansion around the saturation density po'. 

i?(p,0) = -av + -f%(p-po)' + --- (18) 
while the quadratic symmetry energy S2{p) is expanded as: 

S2ip) ^ a4 + ^{p - Po) + ^^{p- Pof -\ (19) 

Po l°Po 

Two key parameters are the symmetry energy at saturation, 04, and the linear 
density dependence, po- (Note: 5*4 (p) is unimportant numerically.) 
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Fig. 10. Skin thickness in ^"^Pb versus the equilibrium symmetry energy 



A study of a wide range of mean- field models, both nonrelativistic Skyrme 
and relativistic mean- field (meson-based and point-coupling) , show strong cor- 
relations between the skin thickness in ^"^Pb and the symmetry energy pa- 
rameters. In Figs. llOland llll this correlation is shown explicitly. While there 
is a large spread in the values of S2{po) = ft4j if we plot the spread (i.e., the 
standard deviation) as a function of p, we see that the spread is smallest at 
a density corresponding to the average density in heavy nuclei (see Fig. I12|l . 
This result is consistent with our conclusion from power counting that only 
one parameter is determined: it is the symmetry energy at some lower density. 

Since pq is a pressure that pushes neutrons out against the surface tension 
in ^'^^Pb, it is instrumental in determining the skin thickness How well do 
we know its value? A comparison of symmetry energies predicted by realistic 
non-relativistic potentials was made in |^ with the result in Fig.^l Reading 
off the slopes implies that pQ « 2.1 MeV/fm with only a 10% spread. From 
relativistic RBHF calculations, pq « 3.3 MeV/fm~'^. In either case, there is 
a strong discrepancy with conventional meson mean-field models, which have 
Pa > 4. 

This conclusion is supported by Danielewicz |42| . who has analyzed the 
relationship in terms of surface vs. volume symmetry energies. Using mass 
formula constraints in conjunction with skin thickness analysis, he concludes 
that 0.56 < iS'2(po/2)/a4 < 0.81, which allows the models between the vertical 
lines in Fig. 1141 Once again, conventional relativistic models are excluded. 
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Fig. 12. The spread (standard deviation) of S^^p) [see Il7ll 1 for many mean-field 
models fit to nuclear bulk properties |88|. 



There are several ways out. Point coupling |45[ I46| . DDRH j37|0Hl> and 
new relativistic mean-field models (with mesons) 031 get lower values of po 
by adding isovector parameters (implicitly, in the case of DDRH). However, 
they achieve lower po from an unnatural cancellation between different orders 
according to the present power counting |38|. This doesn't mean it is wrong. 
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but it means we haven't captured the physics of the symmetry energy. For 
example, it might be pion physics with a different functional form. These 
results imply that we should revisit the symmetry energy in covariant models 
within a chiral EFT. 
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4.2 Guidance from DFT for Solid-State or Molecular Systems 

The density functional theory (DFT) perspective helps to explain successful 
"mean-field" phenomenology, but DFT existence theorems are not construc- 
tive. NDA power counting gives guidance on truncation but not on the analytic 
form of functionals. One approach toward more accurate energy functionals 
would be to emulate the quantum chemists. 

The Hohenberg-Kohn free-energy for an inhomogeneous electron gas in an 
external potential takes the form (ri(x) is the charge density here) 

FHK[n(x)] = F„i[n(x)] + Yjd^x cPx' |^ _ + ^xc[n(x)] , (20) 

where the direct Coulomb piece is explicit and the Kohn-Sham non-interacting 
functional is evaluated using 

A N 
i^ni[n(x)] = j d^a;n(x)weff(x) . (21) 

The '4>iS are solutions to 

(-^V^ +i;eff(x)j^,(x) = e,V,(x) (22) 
with the property that the exact density of the full system is given by 

N 

n(x) =^|^,(x)|2 , (23) 

1=1 

which is the density of a noninteracting system with effective potential 

Vcs ^ V - e4> + Uxc , where -yxc(x) = 6E^c/Sn{-x.) . (24) 

V is the external potential, and (/) is the Hartree potential. 

Solving for the orbitals and energies is the easy part. The hard part is 
finding a reliable approximation to the exchange-correlation functional E^c- 
The approach used by quantum chemists uses the local density approximation 

£;xc[n(x)] « / d^xSM^)) , (25) 



where the energy density fxc is fit to a Monte Carlo calculation of the homo- 
geneous electron gas (see Fig. I15f) . In practice, one uses parametric formulas, 
such as 

£^^{n)/n = -0.458/r^ - 0.0666G(r^/11.4) (26) 

with 
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Fig. 15. The exchange-correlation functional for a uniform system (r^ oc l/fep) |51| . 



G{x) = U{l + xf log(l + x-i) -x^ + ]^x-^- 



(27) 



(the ones in common use are more sophisticated). The calculational procedure, 
which includes correlations through E^c , has the form of a "naive" mean- field 
approach with the additional potential 



dn 



(28) 



n— n(x) 



An example of a successful atomic LDA calculation is shown in Fig. E| Sig- 
nificant improvements to the LDA are made using (generalized) gradient ex- 
pansions |49[ I50| . The DDRH approach to relativistic energy functionals for 
nuclei, in which the density dependence of a nuclear matter calculation is 
incorporated through density-dependent couplings |47l HHj , is in the spirit of 
these LDA plus gradient expansion approaches. 




Fig. 16. The charge density for argon atom calculated using DFT (in the LDA) 
and the Thomas- Fermi approximation |51| . 
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While there are many successes to be found, it is also clear that the situ- 
ation is not entirely satisfactory. From A Chemist's Guide to DFT |52| : 

"To many, the success of DFT appeared somewhat miraculous, and 
maybe even unjust and unjustified. Unjust in view of the easy achieve- 
ment of accuracy that was so hard to come by in the wave function 
based methods. And unjustified it appeared to those who doubted the 
soundness of the theoretical foundations." 

From Argaman and Makov's article 51 (highly recommended!): 

"It is important to stress that all practical applications of DFT rest 
on essentially uncontrolled approximations, such as the LDA ..." 

and from the authors of the state-of-the-art gradient approximations |49| : 

"Some say that 'there is no systematic way to construct density func- 
tional approximations.' But there are more or less systematic ways, 
and the approach taken . . . here is one of the former." 

Thus we are motivated to consider a more systematic approach: EFT! 



4.3 DFT in Effective Field Theory Form 

It is particularly natural to consider DFT in thermodynamic terms, in anal- 
ogy to finding via Legendre transformation the energy as a function of particle 
number from the thermodynamic potential as a function of the chemical po- 
tential This analogy implies the effective action formalism will be well 
suited. The EFT will provide power counting for diagrams and gradient expan- 
sions, thereby systematically approximating a model-independent functional. 
The EFT tells us how to renormalizes divergence consistently and gives insight 
into analytic structure of the functional. 

As a laboratory for exploring DFT/EFT, we can consider a dilute "nat- 
ural" Fermi gas in a confining potential, which is realized experimentally by 
cold atomic Fermi gases in optical traps. Natural in this context means that 
the scattering length is comparable to the range of the potential (rather than 
being fine tuned to a large value) . The EFT Lagrangian takes the form 

^^f' = ^4t + ^ - T^^'^)' - T^^'^)' + ■ ■ ■ ('^^ 

where the parameters can be determined from free-space scattering (or by 
finite density fits); dimensional regularization with minimal subtraction is 
particularly advantageous. The EFT for the uniform system (no external po- 
tential) predicts the energy (and other observables) simply and efficiently 
through a diagrammatic expansion, with a power counting that tells us what 
diagram to calculate at each order in the expansion. [See 53 for details.] The 
first terms in the diagrammatic expansion are (the solid dot is Co = 47ras/M, 
with fls the scattering length): 
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which yields the expansion for the energy density (with degeneracy g): 



2 



C 

c = n 

2M 



l^^^' l){|^(fcFa.) + ^(11 - 21n2)(fcpa,)2 \ + 



(30) 



in powers of kpUs- (Note this is not a power series; see |HS1-) The point here 
is that we have a systematic expansion: LO, NLO, NNLO, and so on, just as 
in the chiral NN case. To put the gas in a harmonic or gaussian trap with 
external potential i'(x), we simply add a v{x.)p{x) term to the Lagrangian. 

The effective action formalism starts with a path integral for the generating 
functional with an additional source J coupled to the density: 

Z[J] = e^'^W = / V^'fViP e*J<^"^[^-J(^)'/'t(^)'/'(^)] (31) 



(Different external sources can be used.) The density is found by a functional 
derivative with respect to J: 

pix) ^ = ■ (32) 

oJ(x) 

The crucial step is the inversion of H32|l: that is, to find J[p]. Given that, the 
effective action r[p] follows from a functional Legendre transformation: 

r[p] = W[J] - J d*x J{x)p{x) . (33) 

The functional r[p] is proportional (with a trivial time factor) to the ground- 
state energy functional E[p] that we seek. The density p(x) is determined by 
stationarity: 

5r \p\ 

= —J(x) =^ in the groimd state. (34) 

dp{x) 

That's density functional theory in a nutshell! 

Two approaches to invert the Legendre transformation have been pro- 
posed. They correspond to generating point coupling or meson-based models. 
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We'll outline the first approach, which is the inversion method of Fukuda et 
al. It is a generalization of the Kohn-Luttinger inversion procedure 

|56|. The idea is to use an order-by-order matching in a counting parameter A, 
which we identify with an EFT expansion parameter (e.g., 1/A^). We expand 
each of the ingredients of H34(l : 



J[p, X] = Jo[p] 4 
W[J,X] = Wo[J] 

r[p,x] = ro[p]^ 



xMp] + ^^Mp] + -- 

^ XWi[J] + X^W2[J] 



(35) 
(36) 
(37) 



and equate each order in A. Zeroth order corresponds to a noninteracting 
system with potential Jq{x): 



ro[p]^Wo[Jo]- Jd''xJo{x)p{a 



p{x) = 



dJo{x) 



(38) 



This is the Kohn-Sham system with the exact density! That is, Jo is the poten- 
tial that reproduces the exact interacting density in a non-interacting system. 
We diagonalize Wo[Jo] by introducing Kohn-Sham orbitals, which yields a 
sum of eigenvalues. Finally, we find Jg for the ground state by completing a 
self-consistency loop: 



Jn 



Wo 



Jo{x) = 



^ 6p{x) 



(39) 



We can contrast the conventional expansion of the propagator: 



3 



+ 




S*(x,x';w) 



which involves a non-local, state dependent i7*(x, x'; cj), with the local Jo(x): 



Jo(x) 




(perms.) 



(perms.) + 



To carry out this construction there are new Feynman rules that involve this 
"inverse density-density correlator," which is indicated by a double line: 

r... - cx3+(A}+oCX>+oOOo + -- 
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Fig. 17. Density profiles for a dilute Fermi gas in a harmonic trap for successive 
approximations. 



Complete details are available in j57j . 

If we restrict the application to closed shells for simplicity, we just specify 
the spin degeneracy g and the Fermi harmonic oscillator shell (the last occu- 
pied shell) Np, which then determines the number A of atoms. The iteration 
procedure for the dilute Fermi gas in a harmonic trap is as follows: 

1. guess an initial density profile p{r) (e.g., Thomas- Fermi) , 

2. evaluate the local single-particle potential Vcs{t) = v(r) — Jo{r), 

3. solve for the lowest A states (including degeneracies) {tpa,ia}' 

^ '"''^''''^'^ "^"^^^ " e«V'a(x) , (40) 

4. compute a new density p{r) — IV'ct(^)Pi other observables are func- 
tional of {ipa^^a}, 

5. repeat 2.-4. until changes are small ("self-consistent"). 

This procedure is simple to implement as a computer program and will be 
familiar to anyone who has performed a relativistic mean-field calculation for 
finite nuclei. 

Some representative results are given in Figs. 1171 andllSI Note the conver- 
gence of the densities and energies as we increase the order of the calculation. 
In Fig. El we show the contribution to the energy per particle at LO, NLO, 
and NNLO for three different systems. We estimate the expected contribu- 
tion by NDA from the Hartree terms in the functional (which appear with 
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Fig. 18. Density profiles for a dilute Fermi gas in a harmonic trap for successive 
approximations. 



simple powers of the density) and then scale correction terms by the expan- 
sion parameter {kF)as. The power counting still works! This is an excellent 
laboratory to explore such issues. For example, the rightmost circle shows a 
large discrepancy between predicted and actual values. A closer examination 
reveals an "accidental" cancellation (that goes away for different parameter 
choices!). 
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Fig. 19. Energy estimates (squares) and actual values (circles) for a dilute system. 
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Although we did not discuss covariant DFT here, Schmid et al. have al- 
ready laid the foundation [SS] ■ The generalization of DFT to other sources and 
time dependence is straightforward (cf. spin-density functionals for Coulomb 
systems). For the covariant case, one treats i/'iV'j ^ matrix and introduces 
a source jij coupled to it, which yields a functional of the scalar density and 
vector current, as described in Sect. 3.3. 

5 Building the Next Generation Theories 

The successes of chiral EFT for nucleon-nucleon scattering and the few-body 
problem have raised the bar for models of nuclear structure. For the next gen- 
eration, we must justify what we do and not just mumble "effective theory" to 
account for phcnomcnological successes. This requires that we develop and test 
systematic power counting schemes and reexamine and sharpen arguments for 
a covariant approach. To meet the challenge, the next generation of covari- 
ant models must move away from model dependence. This implies turning 
away from the "minimal model" philosophy and including all terms allowed 
by symmetries (up to redundancies). It means requiring error estimates and 
a controlled expansion in density, asymmetry, and momentum transfer. And 
it means that the results should exhibit regulator (e.g., cutoff) independence 
up to expected truncation errors. 

There has been significant progress made toward building the next gen- 
eration of covariant approaches to nuclear structure (e.g., see B.D. Serot's 
lecture). But there is a long list of loose ends to address. Some of the physics 
goals and open issues are 

• a reliable description of nuclei far from stability; 

• controlled equation of state and transport properties of neutron stars; 

• unified description of collective excitations; 

• currents, both elastic and transition currents to excited states; 

• improved treatment of isovector physics; 

• consistent incorporation of chiral physics; 

• consistent treatment of pairing in covariant density functionals; 

• establishing more direct connections to QCD (e.g., quark mass dependence 
of equilibrium properties). 

The framework we have proposed to address these goals and issues is a merger 
of effective field theory and density functional theory. The idea is to embed 
the successful relativistic mean-field phenomenology into this framework, as 
traditional NN interactions are taken over by chiral EFT. The basic structure 
and successes are preserved, but model dependence is removed. 

There are many challenges. For the covariant many-body EFT, one needs 
a consistent application of infrared regularization (or equivalent) applied at 
finite density. This is likely to mean the use of cutoff regularization. A pos- 
sible laboratory for these issues is a dilute Fermi gas with A Afjv, where 
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in a finite system the spin-orbit force can be directly explored. In addition, 
chiral EFT expansions must bo translated to covariant form. For the covariant 
DFT, one needs systematic derivative expansions (analogous to past density 
matrix expansions), the incorporation of long-range forces and pairing, and 
a consistent extension to timc-dcpcndcnt DFT. Most importantly, we need a 
more complete power counting scheme for energy functionals. 

There is growing hope for an ultimate path to ab initio calculations of nu- 
clei, starting from quantum chromodynamics (QCD). Nonperturbativc lattice 
QCD could calculate low-energy constants needed for the chiral EFT of NN 
and few-body systems. This EFT in turn determines the analytic structure 
and low-energy constants of a (covariant) many-body EFT, which takes the 
form of a DFT energy functional, applicable across the periodic table. There 
are many gaps, but this dream has become significantly more plausible at all 
levels in just the last few years. 
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